In [1]:
# some plotting instructions
%pylab inline
rcParams['image.aspect'] = 'auto'
style.use('bmh')
In [2]:
from iuvs import io
%autocall 1
In [3]:
import os
files = !ls ~/data/iuvs/level1b/*.gz
for file in files:
print(os.path.basename(file))
In [4]:
l1b = io.L1BReader(files[-1])
The DarkIntegration data shows that there are 2 observations of darks, approx 22 minutes apart:
In [5]:
l1b.DarkIntegration
Out[5]:
Let's look at the original darks. There are 2 available:
In [6]:
l1b.detector_dark.shape
Out[6]:
In [7]:
def compare_darks(dark1, dark2):
fig, ax = subplots(nrows=2, figsize=(10,8))
ax[0].imshow(dark1)
ax[1].imshow(dark2)
In [8]:
compare_darks(*l1b.detector_dark) # this trick puts the first axis of a cube into a function
As one can see, they are strongly different. The last dark looks similar to a usual light image.
In [9]:
l1b.DarkEngineering
Out[9]:
In [9]: